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The vacuum fluctuation (VF) effects on asymmetric nuclear matter are investigated. Masses of 
nucleons and mesons are modified in the nuclear medium by calculating the loop-diagram corrections 
and the density dependence of hadron masses is obtained. The relativistic Lagrangian density with 
the isovector scalar S meson is used to calculate the nuclear equation of state (EOS) in the framework 
of the relativistic mean-field (RMF) approach, the effects of the in-medium hadron masses on the 
properties of neutron stars are finally studied. With the inclusion of the VF corrections, the nuclear 
EOS becomes softer and the neutron star masses are reduced. 
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I. Introduction 



The relativistic mean- field theory (RMF) [l|, 0| , which is one of the main applications of quantum 
hadrodynamics (QHD), is successful in nuclear structure studies 0,0,0]. The nonlinear (NL) Walecka 
model (NLWM), based on the RMF approach, has been extensively used to study the properties of 
nuclear and neutron matter, /3-stable nuclei, and then extended to the drip-line regions 0,0,0,0,0, 
Ell in recent years some authors [ill EE E3, U, 15 1 have stressed the importance of including the 



scalar isovector virtual 5(ao(980)) field in hadronic effective field theories when asymmetric nuclear 
matter is studied. The inclusion of the S meson leads to the structure of relativistic interactions, 
where a balance between an attractive (scalar) and a repulsive (vector) potential exists. The 5 meson 
plays the role in the isospin channel and mainly affects the behavior of the system at high density 
regions and so is of great interest in nuclear astrophysics. 

The properties of nuclear matter at high density play a crucial role for building models of neutron 
stars (NS). Neutron stars are objects of highly compressed matter. The structure of a compact star 
is characterized by its mass and radius, which are obtained from appropriate equation of state (EOS) 
at high densities. The EOS can be derived either from relativistic or potential models. 

In order to describe the medium dependence of nuclear interactions, a density dependent relativistic 
hadron (DDRH) field theory has been recently proposed 16 , l?l Ea Eil ■ Recently the authors 



used the density dependent coupling models with the S meson being included to study the neutron 
stars. They found that the introduction of the <5 meson in the constant coupling model leads to heavier 
neutron stars in a nucleon-lepton picture. The neutron star masses in the density dependent models 
can be reduced when the 6 meson is taken into account. 

The in-medium modification to the masses of er, u>, and p mesons has been studied in experiments 
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and theoretical approaches for a decade 



2Ql, L2JJ, l22l 123J, l24|, Q- Recently the authors of Ref. 
investigated the effect of in-medium meson masses on the properties of the nuclear matter in the 
Walecka model and the effective masses of a and to mesons in the nuclear medium were calculated by 
taking into account the effects of the vacuum fluctuation (VF). In this work we want to see the VF 
effects on asymmetric matter and neutron stars. We also want to clarify the density dependence of 
in-medium nucleon and meson masses. The VF effects are naturally introduced by considering loop 
corrections to the self-energies of in-medium nucleons and mesons. The effective masses of nucleons 
and mesons (cr, u>, p, and 8) in the nuclear medium will be calculated in the VF-RMF model. The 
VF effects on asymmetric matter and neutron stars will be studied. 

This paper is organized as follows. In Sec. II, we derive the in-medium effective masses of nucleons 
and mesons and the EOS for nuclear matter in VF-RMF model. In Sec. Ill, we compare our results 
with those of the NL-RMF model. In Sec. IV, a brief summary is presented. 



II. Hadron effective masses and EOS of nuclear matter 

The relativistic Lagrangian density of the interacting many-particle system consisting of nucleons, 
isoscalar (scalar cr, vector lu), and isovector (scalar S, vector p) mesons used in this work is 

£ = $ i^dn ~ 5^7 M ^ - 5p7 m ^ • b^ - (M - g a (f> - g s r ■ S) ip 



where <j>, uj^, b^, and S represent a, u>, p, and 5 meson fields, respectively, U(4>) = -|ac6 3 + jbcj) 4 is the 
nonlinear potential of the a meson, F„ v = d^oj^ — d v oj^, and G ^ v = d^b v — d v b^. In order to remove 
the divergence in the loop calculation, the counterterm for the Lagrangian density, SC, included in 
Eq. (I) reads 



+y(^) 2 + |(c^) 2 , (2) 



where the parameters cki, /?, and Q (*=1, 2, 3, 4; j=a, u>, p, 6) are fixed by the renormalization 
methods suggested in Refs. 0, Sj- 

The field equations in the RMF approximation are 
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where and p| (the index (subscript or superscript) i 



denotes proton or neutron throughout this paper) are the nucleon and scalar densities, respectively. 
The nucleon and the scalar densities are given by, respectively, 



and 



Pi = -* 



d 4 fc 



TrG 4 (fc), 
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where k Fi is the Fermi momentum of the nucleon and G l (k) is the nucleon propagator in the VF-RMF 
model: 



G\k) = (7^ + M*) 



'U 2 -Mf + ir] + E F . 
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where E F . 



k F . + M* 2 , M* is the nucleon effective masses, and r] is infinitesimal. Here G F (k) 
describes the free propagation of positive and negative energy quasi- nucleons. G l D (k) describes quasi- 
nucleon 'holes' inside the Fermi sea and corrects the propagation of positive energy quasi-nucleons by 
the Pauli exclusion principle. 

The nucleon effective mass without the 8 field in the RMF approach is M* — M — g a <j). When 
the VF effects are considered, the loop-diagram corrections to the self-energy of nucleons and mesons 
shown in Fig. f arc naturally introduced. The nucleon effective mass without the 5 field in the RMF 
approach including the vacuum fluctuations can be calculated from Fig. 1(a). Thus, we have 
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FIG. 1: Loop-diagram corrections to the self-energy of nucleons (a) and mesons (b) in nuclear medium, where 
N denotes nucleon and k is the four momentum of the meson. 

where the second term in the right hand side of the first line in Eq. (7) is given by Fig. 1(a), and the 
third and fourth terms are the contributions from the nonlinear potential of the u meson. 

As well known, the S meson leads to a definite splitting of proton and neutron effective masses, the 
nucleon effective masses with the S meson in the RMF approach are given by, respectively, 



M* p = M- g<r <f>-g s 5 3 , (8) 

and 

M* = M-g a <f> + g s 5 3 . (9) 

The nucleon effective masses with the 5 meson in the RMF approach including the vacuum fluctu- 
ations can be calculated from Fig. 1(a), 
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Thus, we get after the momentum integral 
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where to* (j = c, w, p, J) are the in-medium meson masses. 

The effective mass of the meson (or in-medium meson mass) is defined as the pole position of the 
meson propagator at zero three-momentum (on-shell) or zero four-momentum (off-shell) . We calculate 



the in-medium meson masses in the random phase approximation (RPA) 



27 



. The corresponding 



diagrams are given in Fig. 1(b). Thus, we obtain the in- medium masses of scalar mesons: 
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The expressions of IX, (on-shell and off-shell) are, respectively, 
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The effective masses of vector mesons can be obtained from Fig. 1(b) 



m? = m 2 j +IL jT (<f) (j=u,p), 
where IT^t is the transverse part of the polarization tensor as follows: 
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So the expressions of II^t (on-shell and off-shell) are, respectively, 
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We note that the VF contributions are included in the second summations of ITj(g M ) and IljT'(g^) 
(Eqs. (16), (17) and (20)), the VF contribution in Eq. (21) equals zero. 

The in-medium meson propagator is significantly modified by the interaction of the mesons with 
the nucleons. Clearly, this modification of the meson propagators will change the nucleon effective 
mass as well as the energy density. The meson propagators in the tadpole diagram are calculated at 
zero four momentum transfer. So we must replace the meson mass appearing in the nucleon effective 
mass and the energy density by the meson effective mass defined as 271 ] : 
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Therefore, the energy-momentum tensor in the VF-RMF model can be expressed as 
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Thus the EOS for nuclear matter in the VF-RMF model can be obtained. The energy density is 
given by 
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and the pressure is 
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The density dependence of the nuclear symmetry energy, E sym , is one of the basic properties 
of asymmetric nuclear matter for studying the structure of neutron stars. Empirically, we have 
information on E sym only at the saturation point, where it ranges from 28 to 35 MeV according to 



the nuclear mass table 29(. The nuclear symmetry energy is defined through the expansion of the 
binding energy in terms of the asymmetry parameter a 12j : 
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(27) 



where the binding energy density is defined as E/A = e/p — M, and the asymmetry parameter 

a = (Pn - Pp)/P- 

The nuclear symmetry energy is defined by 
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According to the definition, an explicit expression for the symmetry energy in the VF-RMF model 
is obtained as 



I A. 4- 



114 

4 » r*2 



M* z p 



P 9l M 



2 m s Ef(l + -%C - AAjD) 

Jf V mi 7T Blj ' 

D 



2ir 2 mf Ef l + J^ C _l_J^ D 



p at m 



*2 



4tt 2 mf Ef (1 + Jl c _J_JI D) 2 



12M* 2 ln 



M* 



-7M 2 + 26M(M - M*) - 25(M - M*) 1 



7M 



*2 



(29) 



where 
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As discussed in Refs. 
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nuclear EOS. The incompressibility of nuclear matter is defined by 



the incompressibility K is one of the important ingredients for the 
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where fj, = (e + P)/p is the baryon chemical potential and po is the saturation density. 
We can easily obtain the incompressibility from the EOS in the VF-RMF model as 
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where 
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The final outcome of a supernova explosion can be a neutron star or a black hole. The neutron star 
is believed to evolve from an initially hot protoneutron star. The matter in cold neutron stars is in 
the ground state in nuclear equilibrium. Matter in equilibrium concerning weak interactions is called 
as /3-equilibrium matter. The composition of /3-equilibrium system is determined by the request of 



charge neutrality and chemical-potential equilibrium 3CJ, [32[ . The balance processes for /3-equilibrium 
(npe) system are the following weak reactions: 

n — > p + e~ + P e , (36) 
p + e~ — > n + v e . (37) 

The chemical-potential equilibrium condition for (npe) system can be written as 

p e = Pn - Pp, (38) 



where the electron chemical-potential u e = \/k F + m 2 , kp^ is the electron momentum at the fermion 
level and m e is the electron mass. 



The charge neutrality condition is 



p p = X vPl (39) 



where p e is the electron density, and the proton fraction X p = Z / A = p v j p. 

In the ultra-relativistic limit for noninteracting electrons, the electron density can be expressed as 
a function of its chemical potential 
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Then, we can obtain the relation between the proton fraction X p and the nuclear symmetry energy 
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The EOS for /3-equilibrium (npe) matter can be estimated by using the values of X p , which can be 
obtained by solving Eq. (41). The properties of the neutron stars can be finally studied by solving 
Tolmann-Oppenheimer-Volkov (TOV) equations [33j with the derived nuclear EOS as an input. 



III. Results and discussions 

In order to make a comparison with the NL-RMF model 

00 

, the same saturation properties 

of nuclear matter and hadron masses are listed in Table I, which are used to determine the model 
parameters. The obtained model parameters are presented in Table II with the NL-RMF model 
parameters to geth er for a comparison. The coupling constants are defined as fj = gj/rnj (j = a, u, 
p, 5) in Refs. The parameters of self-interacting terms in Table II are defined as A = ajg\ 

and B = bj g\. 

Table I. Saturation properties of nuclear matter and hadron masses. 



saturation properties 


[12] 


po (/to" 3 ) 


0.16 


E/A (MeV) 


-16.0 


K {MeV) 


240.0 


E sym {MeV) 


31.3 


M*/M 


0.75 


M {MeV) 


939 


m a {MeV) 


550 


to w (MeV) 


783 


to p (MeV) 


770 


m s (MeV) 


980 



Table II. Model Parameters in the VF-RMF and NL-RMF models. 



Parameter 


VF - RMF model 


NL - RMF model [12] 


VFp 


VFp6 


NLp 


NLp5 


9<y 


12.33 


12.33 


8.96 


8.96 




10.52 


10.52 


9.24 


9.24 


9p 


4.01 


6.80 


3.80 


6.93 


gs 


0.00 


7.85 


0.00 


7.85 


A (/to- 1 ) 


0.048 


0.048 


0.033 


0.033 


B 


-0.021 


-0.021 


-0.0048 


-0.0048 



We use the obtained parameters in Table II to complete the calculations of self-consistence in the 
present work. The masses of hadrons (nucleons and mesons) in the medium can be obtained in the 
relativistic mean field approach (RMFA) with the VF effects by calculating the loop-diagrams in Fig. 
1. We first come to the self-consistent calculations of in-medium meson masses. The obtained in- 
medium meson masses are presented in Fig. 2. It is obviously seen from Eqs. (16), (17) and (20), (21) 
that the in-medium mesom masses are related to the asymmetry parameter a. The results in Fig. 2 
are for both cases of a=0.0 (symmetric matter) given by VFp and a=1.0 (asymmetric matter) given 
by VFpS models, respectively. 





Fig. 2(a) shows the in-medium meson effective masses (on-shell: q = 0, c/° 



m~ 



j = a, uj, p, 5) as 



a function of the baryon density in the VF-RMF model. The in-medium modification to the masses 
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21 



22j . The in- medium 



of cr, to, and p mesons have been studied in other theoretical models 
effective mass decrease at the normal density is ^18% for p and u) mesons in the model based on QCD 
sum rule 2l|. The mass decrease is ~20% for u> and p mesons at the normal density according to 



Brown- Rho (BR) scaling law 20j. In our model, the decreases of in- medium meson effective masses 
at the normal density are ~25% for cr, ^20% for lj, and only ^5% for p mesons in the symmetric YFp 
case, respectively. In the VFp5 case, the decreases of p and 8 mesons are 11~13% and 25~27% at 
the normal density, respectively. Most experiments and theoretical approaches indicated a decrease of 
the in-medium meson effective masses around the normal density comparing with the masses at zero 
density 



20. 



21 



22, 



23. 



24. 



25[ . However, in the latest experiment 26[, no significant mass shift for 



the p meson with momenta ranging from 0.8 to 3.0 GeV was observed. Up to now, the experimental 
results have not yet converged, and more work is needed to obtain consistent understanding of the 
in- medium behavior of vector mesons I34j . In general, the medium modifications to the masses of 



mesons are momentum dependent [35| . We note that the on-shell meson effective masses are obtained 
for the mesons at rest in our model. This is not in the momentum range of the CLAS experiment 261 ]. 
Until now, there is no experimental measurement about the in-medium modification to the mass of the 
S meson. Our model indicates a significant decrease of the 6 meson effective mass around the normal 
density. We note that the effective meson masses become to increase at high density regions in the 



VF-RMF model. Unfortunately, the high density regions are beyond the reach of current experiments. 
It will be very interesting to test our prediction in the future experiments. 

Fig. 2(b) shows the in-medium meson effective masses (off-shell: — 0), which are used for the 
calculations of the nuclear EOS, as a function of the baryon density in the VF-RMF model. The 
off-shell meson masses are different from the on-shell meson masses because the four momenta carried 
by the meson propagators in these two situations are different. Since the meson propagators in the 
nucleon self-energy are computed at zero four momentum transfer (see Fig. (Ha)), we have to use the 
off-shell meson masses in the tadpole loop calculation for self-consistency. 
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The nucleon effective masses play an important role in the calculations of the EOS (see Eqs. (25) 
and (26)). We calculate the loop-diagram corrections to the self-energies of nucleons in medium. The 
nucleon effective masses without the S meson can be calculated from Eq. (7). One can see from Eqs. 



(12) and (13) that the presence of the S meson leads to proton and neutron effective mass splitting. 
We present the baryon density dependence of proton and neutron effective masses for different proton 
fractions in the two models for a comparison in Fig. 3. The solid lines in Fig. 3 are the nucleon 
effective mass for symmetric matter (X p =0.5). 

Fig. 3(a) shows that the proton and neutron effective masses given by the VFpS model decrease 
slowly with the increase of the baryon density, at variance with the NLpS model that presents a much 
faster decrease (Fig. 3(b)). This main difference between the VFpS and the NLp<5 models actually 
comes from the in- medium meson masses (see Eqs. (22) and (23)). The in- medium meson masses 
(off-shell) increase with the increase of the baryon density in the VF-RMF model (see Fig. 2(b)). So 
it is natural that the VF effects lead to a slow decrease of the nucleon effective masses as the baryon 
density increases. 
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FIG. 4: The symmetry energy as a function of the baryon density in different models. The inset is the 
corresponding proton fraction. 



The density dependence of the symmetry energy for the VF-RMF and NL-RMF models is presented 
in Fig. 4. We see a similar behavior of E sym at saturation density for both the two models. With the 
increase of the baryon density, the symmetry energy given by the VF-RMF model increases slowly 
comparing with that given by the NL-RMF model. From Fig. 4 we see that the symmetry energy 
with the 5 meson is stiffer than that without the S meson for both the VF-RMF and the NL-RMF 
cases. The symmetry energy in the VFp<5 case is softer than that in the NhpS case. This is due to 
the VF effects. The presence of the S meson affects the symmetry energy and consequently the EOS 
of asymmetric nuclear matter. 

The /3-equilibrium matter is relevant to the composition of the neutron stars. The EOS, pressure 
vs density, for (npe) matter in the VF-RMF and the NL-RMF models is presented together in Fig. 5 
for a comparison. We see that the EOS in the VF-RMF is lower. Due to the VF effects, the EOS of 
asymmetric matter becomes softer. 

In the present work, only two pictures for the neutron star composition are considered: pure neutron 
and /^equilibrium matter, i.e. without strangeness bearing baryons and deconfined quarks (see Refs. 



32, 36]). Furthermore, we limit the constituents to be neutrons, protons, and electrons in the latter 
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FIG. 5: Equation of state for (npe) matter. 

case. In fact, in /3-equilibrium matter, nucleons and electrons indeed dominate at low temperature. 

The structure of neutron stars can be calculated by solving TOV equations. The correlation between 
the neutron star mass and the corresponding radius for the pure neutron and the /3-equilibrium (npe) 
matter by the VF-RMF model are shown in Fig. 6. The obtained maxium mass, corresponding radius 
and central density are reported in Table III. We see from Fig. 6 and Table III that the VF-RMF 
model leads to the decrease of the neutron star masses for both the pure neutron and the (npe) 
matter. However, the NL-RMF model leads to heavier neutron stars (see Refs. [3, This is 

mainly because the EOS of asymmetric matter becomes softer since the symmetry energy becomes 
softer in the VF-RMF model. 
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FIG. 6: The mass of the neutron star as a function of the radius of the neutron star. 



We note that the coupling constant fj = gj/p j in the NL-RMF model is a constant fixed by the 
saturation properties of nuclear matter [ill [l2 |. In order to make a comparison, we define Fj — 
9j/ m j 2 (p) =(J j ^, P, S) in our model, here m* is the off-shell meson mass. We are interested in the 
VF effects on the coupling constants. It is well known that the coupling constants in the DDRH model 
are density dependent (see 16[), whereas the meson masses are constant. In the present model, the 
coupling constants are constant, but the meson masses are density dependent. In order to distinguish, 
we define Fj — g* 2 (p)/m 2 j (j — a, ui, p, g* is the coupling constant) for the DDRHp case. We present 
a comparison between F* and Fj in Fig. 7. We see from Fig. 7 that F* decreases with the increase 



Table III. The maximum mass, the corresponding radius and the central density of the neutron star in 

the VF-RMF model. 





Model 


VF- 


RMF 


neutron star 


properties 


VFp 


VFpS 


pure neutron 


M,s/M Q 


1.82 


2.07 




R(km) 


11.57 


12.49 




pel pa 


6.41 


5.48 


(npe) matter 


Ms/Mq 


1.45 


1.51 




R(km) 


10.22 


10.77 




Pc/P0 


8.98 


8.20 



of the baryon density for both a and u) mesons, which are due to the increase of the off-shell meson 
masses in high density regions. It shows that the variance trends of F* and Fj are roughly the same 
for (7, u>, and p mesons. If we attribute the density dependence of the coupling constants in the DDRH 
model to the VF effects in our model, the density dependence of F* and Fj should have the same 
physical origin. 
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FIG. 7: A comparison between the ratio Fj and Fj. 
IV. Summary 

The VF corrections are investigated in the framework of the RMF approximation by using the 
relativistic Lagrangian density with the S field in this work. By taking into account the loop corrections 
to the self-energies of the in-medium nucleons and mesons, the VF effects are naturally introduced 
into the RMF model. In order to make a comparison with the NL-RMF model, the same saturation 
properties of nuclear matter are used to determine the parameters of the VF model. We calculate 
the contributions from self-energy diagrams to the masses of the nucleons and mesons. The effective 
masses of the nucleons and mesons, especially the p and 8 mesons, in the nuclear medium are obtained. 
We find that the nucleon effective masses decrease more slowly with the increase of the baryon density 




comparing with the NL-RMF case. We also find that the dependence of the off-shell meson masses 
on the density in the medium is different from that of the on shell meson masses (see Fig. 2). The 
effects of the in-medium hadron masses on the nuclear EOS and the properties of the neutron stars 
are studied. The VF effects lead to softness of the symmetry energy and the EOS of asymmetric 
matter. Due to such softness, the neutron star masses are reduced quite a lot. This indicates that the 
VF effects on the neutron stars are important. We see from Fig. 2(b) that the off-shell in- medium 
meson masses increase with the increase of the baryon density. The off-shell in-medium meson masses 
can be used to calculate F*. The variance trends of F* for a, u>, and p mesons are roughly consistent 
with Fj in the DDRHp case 16|. The density dependence behavior of the off-shell in- medium meson 
masses is very interesting indeed. 

In the present work, we work in the Hartree approximation in which we only consider the dominant 
contributions from tadpole diagrams to the nucleon self-energy when we investigate the VF effects. 
In fact, the exchange diagram contributions can only provide small corrections to the EOS in the 
RMF approach at high densities It is well known that the symmetries of the infinite nuclear 
matter system can simplify the mean- field Lagrangian considerably. As pointed out in Ref. [2J , when 
translational and rotational invariance of the nuclear matter is taken into account, the expectation 
values of all three-vector fields must vanish. Therefore, the expectation value of in Eq. (I) is 
zero and nonlinear p meson interactions (three or four p vertices) do not appear. Furthermore, since 
we only take into account the tadpole diagrams, only the neutral iso-vector meson (p°) is involved in 
the nucleon self-energy diagrams even in the asymmetric nuclear matter (p^ could contribute to the 
exchange diagrams which are ignored in the Hartree approximation). 

It is suggested that the tensor coupling to the nucleon should be taken into account for the study 
of the vector mesons in nuclear medium [37]. In the present work, since we focus our attention on 
the VF effects on the properties of neutron stars, we do not include the tensor coupling effects in the 
calculations for the effective masses of the vector mesons. This will be studied in the future work, 
especially for the study of the p meson masses in the medium. Furthermore, the future study of 
the VF effects on the meson- nucleon coupling constants will be carried out. In addition, more careful 
study of the high density behavior of the meson- nucleon effective couplings, especially gs, is important 
and attractive. 
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